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NOMENCLATURE 


C. = heat capacity of jth region 

j 

T = temperature of jth region 

T = temperature of kth region 

K 

(T. = derivative of T. with respect to time) 

J J 

C = over-all conduction coefficient between regions k and j 

R. . = over-all radiant coefficient for net radiative transfer from 

kj . , . 

region k to j 

q. = internal power dissipationin the jth region 


g 

A. = effective area of surface j to sun 
J J 

A. = total at-ea of surface j radiating to space 
€ = infrared emissivity of surface j 

a . = absorptivity of surface j with respect to solar insolation 


S 

a = 



d, b, 

k. = 
J 

(PC P 1 

H. = 
J 


solar constant 

Stefan- Boltzmann constant 

classical geometry factor for net radiation exchange based on 
Lambert's cosine law for diffuse radiation 

r i, r 2 , r 3 , r 4 , L, W, h = particular linear dimensions (Fig. Z) 

thermal conductivity of jth region 

= density-specific heat product for jth region (volumetric heat 
capacity) 

s 

solar radiation and internal generation ( = q,+ a,A .S) 


t = time 


IV 



TECHNICAL MEMORANDUM X-53346 


A SET OF EXPERIMENTS IN, 

THERMAL SIMILITUDE 

SUMMARY 

The analysis and results for two sets of experiments in thermal 
scale modeling are presented. The prototype and model consist of a 
plate, cylinder, and sphere exchanging thermal energy by radiation 
only. They were located relative to one another in an unsymmetrical 
arrangement. 

The experimental results generally confirm the modeling rules, 
with some exceptions in the details, due largely to an a priori assump- 
tion made regarding the volume partitioning of the objects to be modeled. 



INTRODUCTION 


* 


Previously derived scale modeling laws [l] were used as the basis 
for the design of experiments in thermal similitude. The purpose of the 
experiments was to test the theoretical laws for the transient case where 
several bodies with geometrical shapes common to space vehicles ex- 
change energy solely by thermal radiation. The three bodies, a sphere, 
a cylinder, and a flat plate, were in a particular asymmetrical, 
geometrical arrangement (Figure 1). 

Experiments were conducted on both a full-scale prototype and a 
model, since the primary objective was to test the scaling laws. The 
main overall dimensions of the prototype were scaled a priori by j and 
minor dimensions were deliberately distorted to satisfy the scaling re- 
lations. 

An important secondary objective of the experiments was to dis- 
cover practical difficulties in modeling and to develop techniques to 
overcome them. 

Analysis for the comparison of the prototype and model tempera- 
tures was made on an IBM 7090 computer, and the graphs were machine 
plotted using a Stromberg-Carlson SC 4020 plotter. 

Some modeling work done by Adkins that is very similar to this 
has been reported in the literature [ 2] . 

This report includes two sets of experiments which were conducted 
several months apart. 
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DESCRIPTIVE EQUATIONS 


Since there was no conductive or convective heat exchange between 
the objects, only radiative exchange will be considered. They are assum- 
ed to be exposed to external radiation equivalent to solar radiation and to 
have heat dissipated internally. Part of the objective was to use the cru- 
dest possible physical partition of the bodies in order to determine the 
limits on crudeness for modeling. Toward this end, each object is re- 
presented by a single temperature. The difference equations describing 
their transient behavior under these conditions are [3] 



since T 4 represents the chamber wall, which is maintained with liquid 
nitrogen at a fixed temperature. 


SCALING RELATIONS AND MODEL DESIGN 


The general similarity criteria which includes this special case 
is given by Eq. (4) in Reference [l] . From this list, with the additional 

definition 


s 

H. = q. + u.A .S , j = 1 , 2, 3 . 
3 J 3 3 


(2) 


3 



the following independent set of scaling ratios is obtained: 




T. € ctA T 3 t 

_1 _J L 

T, ’ C. 

k J 


R, ,T 3 t 
k J J 
C 

J 


H.t 

C.T. • 
J J 


( 3 ) 


If electrical resistance heaters are used to obtain the simulated 
heating effects of solar insolation, and if the space chamber is regarded 
as incorporated in the third ratio with its appropriate exchange coeffi- 
cient, then the modeling ratios for the problem are 


T. R, .T t 
_1 > k J 1 
T, C. 



(j, k 


1, 


4) 


(4) 


The detailed dimensional notation used is shown in Figure Z. 

All major external dimensions are scaled by \ . It is assumed that the 
materials are not changed from prototype to model, that the radiation 
geometry factors remain unchanged, that the temperatures of the model 
at a particular time are equal to the corresponding prototype tempera- 
tures at the same time, and that the thickness of the cylinder end caps 
are not changed from prototype to model. These a priori requirements 
are 

T. = T. t' = t ( p c p ) .' = ( pc p ). 

J J P J P J 



2 

w" 

b" 


1 

2 v z 

h" = 

£h 

U 

r 4 = 

1 

2 r 4 

= b 




(j = 1,2, 3) 


(5) 
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These requirements, together with the scaling laws [4] , impose 
geometric distortion in the minor dimensions. The results, summarized 
in Table I, were used as the basis for the detailed model design. 

It is economical to fabricate whenever possible from standard 
gage material. However, review of the modeling rules in Table I re- 
veals that the minor dimension distortions (thicknesses) constitute one 
of the main aspects. In addition, standard gage thicknesses have an 
average tolerance of ± 0 . 005 inches, (1. 27 x 10~ 4 meters) and there is 
often non-uniformity within a particular stock sheet. To minimize 
these difficulties, it was decided that the larger objects (prototype) 
would be fabricated from gage stock, but that the thicknesses of the 
stock sheet used would be measured to the smallest possible 
tolerance prior to use. Then, these measured thicknesses would be 
used in the modeling relations to predict the dimensions required for 
the models. The models, which are small, would then be fabricated 
by machining and milling to the tolerance of +0. 001 inch (2. 54 x 10 meters). 
This appeared to be the most economical approach without compromising 
the modeling laws. The resulting dimensions are shown in Table II. 

Electrical resistance heaters were used to represent internal 
dissipation and solar insolation on the plate. The heat capacity of the 
heater would upset the modeling criteria if it is relatively large com- 
pared to the heat capacity of the plate. Devices for attaching the heater 
to the plate also contribute substantially to the heat capacity of the sys- 
tem. However, the heater in all its aspects could also be modeled. It 
was decided that it would be best to represent the plate by two parallel 
sheets, each having equal thickness. The heater would then be mounted 
between the plates by carefully designed connectors, and the heater 
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TABLE I 


Explicit Modeling Criteria 


Quantity 

Model/Prototype 

Remarks 

T. 

J 

i 

a priori determined 

t 

i 

a priori determined 

F kj 

i 

a priori determined 


r 

a priori, material property 

L 

1 

2 

a priori, major dimension 

W 

1 

2 

a priori, major dimension 

r 2 

1 

2 

a priori, major dimension 

r 4 

1 

2 

a priori, major dimension 

h 

1 

2 

% a priori, major dimension 

q j 

1 

4 


d 

1 

distorted (derived) 

b 

1 

distorted (a priori) 

r ! 

1 k 

distorted (derived) 

r 3 

^i k 
2 3 

distorted (derived) 


where 



and 


k 


3 


1 


l 

" 2 



1 

2 
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TABLE II 


Detailed Dimensions of Prototype and Model 


Dimensions 

Prototype 

Model 


Stock 

Measured 

Measured 


Inches 

(Meters) 

Inches (Meters) 

Inches 

(Meters) 

L 



28. 364 (. 72042) 

14. 182 

(. 3602) 

W 



15. 709 (. 3989) 

7. 854 

(. 1995) 

d 

0. 1285 

(. 003264) 

0. 257 (. 006528)* 

0. 1285 

(. 003264) 

b 

0. 064 

(. 001626) 

0. 065 (. 001651) 

0. 065 

(. 001651) 

h 



7. 726 (. 1962) 

3. 832 

(. 09733) 

r 

l 

5. 866 

(. 1490) 

5. 880t(. 1493) 

2. 877 

(. 07307) 

r 

2 

6. 000 

(. 1524) 


3. 000 

(. 07620) 

r 

3 

1. 706 

(. 04333) 

1. 707*(. 04336) 

0. 701 

(. 001780) 

r 

4 

1. 964 

(. 04988) 


0. 982 

(. 02494) 


*r 2 - r x was measured 0. 120 inches (. 003048 m), r 4 - r 3 was 
measured 0. 257 inches (. 006528 m). 

^Measured value was 2d 
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itself would have a very small characteristic time and low heat capacity. 
The sandwich arrangement was to be used for both the prototype and 
model. This arrangement would have allowed an approximate scaling 
of the heaters according to the following reasoning. By the modeling 
rules, the area of one side of the model plate would be 4 that of the pro- 
totype. Assuming that the same heater material in the same wiring 
arrangement is used for both the prototype and model, the heat capacity 
of the heater will be directly proportional to the area of one side of the 
plate in both. So, the heater heat capacity for the model will be ^ of 
that of the prototype. Since the internal generation is also scaled by \ , 
the ratio will not be affected by the heater. In all of the ratios 

R, ,/C , C is scaled by 4, since d* = d ; there is also ai overall scal- 

k 1 1 1 

ing of the numerator, so these ratios would not be affected either. It 
turned out that the sandwich construction was too difficult, especially 
with regard to the heater construction and mounting. Therefore, grooves 
were made in the plates. The wiring was routed through the grooves 
which were then filled with a special thermally conductive and electrically 
insulating cement. Roughly speaking, the same reasoning applies as for 
the intended sandwich construction, except that in making the grooves the 
plate material was removed before the cement material was added. In 
this respect, it was to a great extent fortuitous that the detailed modeling 
criteria presented in Table I was used. 

Welding was used to fabricate the cylinder and to mate hemi- 
spheres. The amount of material in these welds contributes to the heat 
capacity. In order to minimize error from this source, all welds were 
ground to the contour of the external surfaces and care was taken to pre- 
vent excessive accumulation of weld material on the inside. 
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The thickness of thermal coatings could contribute significantly 
to the heat capacity, especially in the very small models. This problem 
was solved by smoking the surfaces of both prototype and model with car- 
bon black produced by an acetylene torch. Resulting emittances were 
0. 98 to 0. 99, according to measurements. 

Another anticipated difficulty was the crudeness of the volume 
partition of the objects to be modeled. The finer the partition, the closer 
the basic equations [l] (and consequently, the modeling ratios) represent 
the objects thermal properties. On the other hand, as the partition be- 
comes finer, the number of ratios that must be satisfied increases rap- 
idly and the model design becomes more complicated. The crudest pos- 
sible partition was used here, and it was anticipated that this would cause 
difficulty in analyzing the experimental data if there were large tempera- 
ture differences around the object. Since the sphere for the first set of 
tests was made of stainless steel, it was the object most likely to exhibit 
this effect. 

EXPERIMENTAL APPARATUS AND PROCEDURES 

Both prototypes and models were suspended in the vacuum cham- 
ber by long, low-conductivity, small-diameter cables to prevent con- 
ductive exchange between the objects and the chamber walls. The rec- 
tangular chamber was 1. 524 x 1. 524 meters in cross-section and 3. 658 
meters in length. Pressure in the chamber during all experiments was 

- 4 

10 mm Hg or less. The walls were covered with cryopanels which pro- 
vided a low temperature inclosure. In addition, both prototype and model 
occupied the same place in the chamber and were oriented in the same 
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way to minimize chamber influence on the test objects. The models as 
they appeared when installed in the chamber are shown in Figure 3. 

Two sets of experiments were conducted. In the first set, stain- 
less steel spheres and a plate with heaters installed were used. In the 
second set, changes were introduced in the sphere to obtain a more uni- 
form temperature; the sphere in which heaters were installed was fab- 
ricated from aluminum alloy. 

The three test objects had a total of 34 thermocouples for the 
first set of tests and 25 for the second set. Seven thermocouples were 
used to monitor the vacuum chamber cryopanel temperatures. The in- 
dividual thermocouple identification numbers and locations are listed in 
Table III and Figure 4. On Table III, thermocouple number 1 for the 
first experiment and number 1 for the second experiment are at the 
same location on the sphere. This is true for thermocouple numbers 
1 through 11, but not for any thereafter . For example, thermocouple 
number 13 for the first experiment and number 12 for the second are at 
the same location on the sphere. The locations can be seen by referring 
to Figure 4 in which the thermocouple numbers refer to the first experi- 
ment only. A minimum number of carefully located thermocouples was 
specified in order to minimize the heat conduction through the wires and 
still obtain adequate information on the temperature differences on op- 
posite sides of the sphere and the cylinder. 

Power loads during the power-on phase of the experiments are 
listed in Table IV. 
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TABLE III 


v 


Individual Thermocouple Identification Numbers and Locations 



Thermocouple 

Number 

Object 

First 

Second 

Experiment 

Experiment 


i 

i 


2 

2 


3 

3 


4 

4 


5 

5 


6 

6 


7 

7 

Sphere 

8 

8 


9 

9 


10 

10 


11 

11 


12 

-- 


13 

12 


14 

13 


15 

14 


16 

15 


17 

16 


18 

17 


19 


Cylinder 

20 



21 



22 



23 

18 


24 

19 

1 

25 

20 


26 

21 


27 

22 


28 

23 

Flat Plate 

29 

24 


30 



31 



32 

25 


33 



34 



Note: Thermocouples 1 through 11 are at the same location in both experi- 
ments. Thereafter, the numbers vary. For example, thermocouple #13 for 
the first experiment and thermocouple #12 for the second experiment are 
in the same location on the sphere, etc. 
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TABLE IV 


Power Loads 
(watts) 


Object 

Prototype 

Model 


Set 1 

Set 2 

Set 1 

Set 2 

Plate 

1053 

1053 

264 

264 

Sphere 

0 

400 

0 

100 

Cylinder 

0 

0 

0 

0 


An important aspect of modeling in the transient case is the 
initial conditions of the experiments. The temperature distribution 
for both model and prototype must be the same at the start of the test 
(or, in general, if temperatures are scaled, then an analogous state- 
ment holds for homologous temperatures at zero time). The difficulties 
arise when considering the duplication of cryopanel cool-down and heater 
turn-on. In order to eliminate possible differences in individual cryo- 
panel cool-down curves due to different ambient starting temperatures, 
nitrogen flow rates, and manifolding between panels, the following ex- 
periment sequence was specified and used. A vacuum was established, 
then for both prototype and model, heaters and nitrogen for the panels 
were activated simultaneously and the systems were allowed to go to 
practical equilibrium. (If the heaters were not also activated, the 

objects under experiment would go slowly to panel temperature, taking 
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much more time and possibly causing difficulty with heaters, paints, 
and instrumentation at the low temperatures that would be obtained). 

A detailed analysis of various conditions is given in Ref. [3] . These 
equilibrium conditions were considered to be the initial conditions for 
the experiments. Then, the actual experiment sequence was cooling 
from these initial conditions (heaters suddenly "off") for two hours, 
then heating (heaters suddenly "on") for two hours. The periodic nature 
of these experiments has the inherent advantage of eliminating or mini- 
mizing the initial conditions which act as transients in the system, and 
which are eventually damped out by the thermal inertia of the system. 


RESULTS AND DISCUSSION 

Figures 5-l6 show the experimental results for the first set of 
experiments. In all such graphs, temperature is plotted as a function 
of time and each prototype and corresponding model temperature is 
plotted on the same graph. The results for all cylinder and plate tem- 
peratures are excellent, and the model temperatures would be entirely 
adequate for predicting prototype temperatures. Two of the tempera- 
tures for the sphere (Figures 6 and 8) are acceptable, the differences 
between prototype and corresponding model temperatures being less 

o 

than 8 K. However, two of the model temperatures for the sphere 
(Figures 5 and 7) are not good predictors of the corresponding proto- 
type temperatures. The reason for this is clear. One of them, thermo- 
couple number 3 (Figure 7), is the one nearest the heated plate, and the 
other, thermocouple number 1 (Figure 5), is the one farthest away 
from the heated plate. The low conductivity of the stainless steel sphere, 

together with the asymmetric heating by radiation from the plate, is not 
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consistent with the crude volume partition of the sphere. In other 
words, the a priori assumption that the model sphere can be designed 
using a single representative temperature is too crude. This is shown 
clearly in Figure 29, in which all sphere temperatures are plotted on 
one graph. It should be noted that all curves have the proper general 
shape. 

The experimental results for the second set of experiments are 
shown in Figures 17-28. Here the spheres were constructed of alumi- 
num and had some power dissipation in them (Table IV) in an attempt 
to obtain more uniform sphere temperatures. The expected improve- 
ment was obtained (compare Figures 29 and 30). The model tempera- 
tures for the sphere may be used to predict the corresponding prototype 
temperatures with good results, as can be seen from Figures 17 through 
20. The cylinder curves (Figures 21-24) are not as good in this experi- 
ment as they were in the first. Again, as can be seen by reference to 
the geometric arrangement (Figures'l and 3), the combined heating of 
the sphere and the plate with asymmetric relative locations makes the 
assumption of a single representative temperature for the cylinder 
somewhat questionable (Figures 31-32) although not nearly as much so 
as for the sphere. In addition, and perhaps more importantly, the ini- 
tial conditions for the prototype and model were not properly established; 
this is shown clearly by the graphs. Furthermore, the results for the 
plate (Figures 25-28) are not nearly as good in the second experiment 
as in the first. This is surprising since both sets of experiments used 
the same prototype and model plates. Furthermore, the reason for 
these apparently anomalous results for plate cannot be due to the fact 
that the sphere was heated in one experiment and not in the other since 

all temperatures over the plate show this same effect. The heating 
14 



portion of the curves could be explained by overheating of the models, 
but this would not explain the cooling portion of the curves. However, 
both heating and cooling curves can be explained by a lower emissivity 
for the model than the prototype. Coatings for these were smoked-on 
carbon black, which handling removes easily, and it is plausible that 
some of the carbon black was inadvertently removed from the model 
plate. The results for the second set of experiments are also shown 
in Tables V, VI, and VII. They are self-explanatory. 


CONCLUSIONS 

1. The modeling laws are generally confirmed by the experi- 
ment. 

2. The volume partitioning used here was too crude; the sphere 
should have been divided into at least four regions. 

3. Electrical resistance heaters are very useful to simulate 
internal energy dissipation but are not entirely satisfactory to simulate 
incident radiation because they necessarily become a part of the objects 
to be modeled, requiring that they also must be modeled. 

4. Considerable care must be used in obtaining proper initial 
conditions if model experiments are to be used to predict prototype 
temperatures in transient cases that are not periodic. 
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TABLE VI 


Fraction of Time in Percent where the Difference Between the Prototype 
and Model Temperatures of the Individual Thermocouple Measurements 
are Equal to or Less than 5, 10, 15, Z0, and Z5 Degrees Kelvin. t 


Thermocouple / ^ 
Number / 

10 

15 

20 

25 

1 

59. 2 

75. 5 

89. 8 

95. 9 

98. 0 

2 

26. 5 

46. 9 

63. 3 

79. 6 

95. 9 

3 

51. 0 

79. 6 

89. 8 

89. 8 

93.9 

4 

55. 1 

73. 5 

89. 8 

93. 9 

93. 9 

5 

44. 9 

57. 1 

69. 4 

79. 6 

89. 8 

6 

51. 0 

65. 3 

85. 7 

89. 8 

91. 8 

7 

51. 0 

75. 5 

83. 7 

85. 7 

87. 8 

8 

55. 1 

71. 4 

85. 7 

89. 8 

89. 8 

10 

53. 1 

67. 3 

87. 8 

89. 8 

91.8 

11 

42. 9 

61. 2 

81.6 

89. 8 

91. 8 

12 

2. 0 

6. 1 

10. 2 

• 14. 3 

14. 3 

13 

0 . 0 

2. 0 

18. 4 

40. 8 

44. 9 

14 

38. 8 

51. 0 

55. 1 

57. 1 

59. 2 

15 

10. 2 

28. 6 

59. 2 

79. 6 

95. 9 

16 

6. 1 

36. 7 

75. 5 

93. 9 

100. 0 

17 

4. 1 

30. 6 

67. 3 

81. 6 

93. 9 

18 

4. 1 

8. 2 

32. 7 

63. 3 

79. 6 

19 

4. 1 

10. 2 

49. 0 

71. 4 

81. 6 

20 

32. 7 

46. 9 

46. 9 

46. 9 

59. 2 

21 

42. 9 

49. 0 

53. 1 

59. 2 

83. 7 

22 

4. 1 

24. 5 

46. 9 

46. 9 

49. 0 

23 

2. 0 

2. 0 

2. 0 

26. 5 

65. 3 

24 

44. 9 

46. 9 

46. 9 

57. 1 

77. 6 

25 

46. 9 

46. 9 

46. 9 

53. 1 

73. 5 


t Determined from discrete points in time corresponding to data 
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Figure 1 - Geometric Arrangement. Numbers shown are based on outside 
dimensions of sphere's radius. Coordinate system shown is 
centered at base of cylinder. 
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(a) FLAT PLATE 


d 




(c) SPHERE 

Figure 2 - Dimensional Notation-Numerical Values Given in Table II 
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igure 3 - Model Installed in Vacuum Chambe 
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Figure 6 - Temperature-Time, Prototype and Model Sphere, Thermocouple 
No, 2 (First Experiment) 
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Figure 7 - Temperature-Time, Prototype and Model Sphere, Thermocouple 
No, 3 (First Experiment) 
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Figure 9 - Temperature-Time, Prototype and Model Cylinder, Thermocouple 
No. 17 (First Experiment) 
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TIME (IN MINUTES) 


Figure 10 - Temperature -Time, Prototype and Model Cylinder, Thermocouple 
No. 19 (First Experiment) 
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Figure 11 - Temperature- Time, Prorotype and Model Cylinder, Thermocouple 
No. 22 (First Experiment) 
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Figure 12 - Temperature-Time, Prototype and Model Cylinder, Thermocouple 

No. 24 (First Experiment) 35 
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Figure 14 - Temperature -Time, Prototype and Model Plate, Thermocouple 
No. 27 (First Experiment) 
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Figure 16 - Temperature-Time, Prototype and Model Plate, Thermocouple 
No. 32 (First Experiment) 
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Figure 17 
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Temperature-Time, Prototype and Model Sphere, Thermocouple 
No. 1 (Second Experiment) 
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Figure 20 - Temperature-Time, Prototype and Model Sphere, Thermocouple 
No. 4 (Second Experiment) 
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Figure 22 - Temperature-Time, Prototype and Model Cylinder, Thermocouple 
No. 17 (Second Experiment) 
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Figure 23 - Temperature-Time, Prototype and Model Cylinder, Thermocouple 
No. 18 (Second Experiment) 
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Figure 24 - Temperature-Time, Prototype and Model Cylinder, Thermocouple 
No. 19 (Second Experiment) 
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Figure 25 - Temperature-Time, Prototype and Model Plate, Thermocouple 
No. 20 (Second Experiment) 
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Figure 27 - Temperature-Time, Prototype and Model Plate, Thermocouple 
No. 24 (Second Experiment) 
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Figure 28 - Temperature-Time, Prototype and Model Plate, Thermocouple 
No. 25 (Second Experiment) 
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Figure 30 - Temperatures for Prototype Sphere for Second Experiment 
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Figure 32 - Temperatures for Prototype Cylinder for Second Experiment 
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